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ABSTRACT 


The  Aeronautics  Department  shock  tube  has  been  developed  and 
instrumented.  The  shock  tube  has  been  used  in  an  experiment  in  which 
the  thermal  conductivity  of  argon  has  been  determined  in  the  temperature 
range  of  1500  -  5000  °K  and  at  relatively  high  pressures,  10  to  30 
atmospheres. 
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CHAPTER  I 

INTRODUCTION 

Although  the  shock  tube  has  been  known  as  an  experimental  device 
in  the  study  of  gas  dynamics  for  the  past  70  years'-  *   since  Paul 
Vieilli  in  France  (1899)  made  first  known  use  of  it,  it  was  not  until 
the  early  1940' s  when  better  methods  and  equipment  were  developed  to 
measure  the  phenomena  occuring  in  the  gases  tested,  that  it  began  to 
enjoy  much  application  or  widespread  use. 

With  the  rapid  technological  development  of  electronic  measuring 
equipment  permitting  the  use  of  fast  response  high  sensitive  instruments 
the  shock  tube  has  become  a  very  useful  experimental  tool  for  study 
and  research  in  the  fields  of  gas  dynamics.  Specific  aDnlications  include 
the  study  of  shock  waves,  shock  fronts  and  interactions,  aerodynamics 
for  the  study  of  the  subsonic,  transonic  and  supersonic  flow  regions, 
as  well  as  a  means  of  driving  hypersonic  wind  tunnels'-  ,  chemistry 
for  the  study  of  high  temperature  gas  reactions'-  J ,  physics  for  study 

of  conversion  of  energy  from  one  form  to  another,  i.e.,  translational 

["21 

to  energy  of  internal  vibration  on  rotation  of  molecules'-  ,  in  the 

field  of  planetary-entry  fluid  physics  connected  with  the  design  of 
the  spacecraft's  heat  shield1-  -1  and  many  other  applications. 

The  purpose  of  the  work  presented  in  this  thesis,  was  first  to 
install  instrumentation  and  leave  in  working  condition  a  newly 
developed  shock  tube  for  the  Department  of  Aeronautics  of  the  Naval 
Postgraduate  School  and  after  that  was  accomplished,  use  it  to 
determine  the  thermal  conductivity  of  a  noble  gas  (argon  in  this 
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case)  in  the  range  from  1500  to  5000  degrees  Kelvin  at  relatively 
high  pressures  i.e.  10  to  50  atmospheres.  Although  several  authors 
have  reported  results  in  this  type  of  experiment"-  ■'''■   ,  they  all 
have  been  done  at  low  pressure  i.e.  around  1  atmosphere  and  data  for 
the  same  experiment  worked  here  was  not  available.  The  value  obtained 
here,  for  the  calibration  parameter  of  the  shock  tube,  3  ,  is  some 
15%  higher  than  similar  parametric  values  verified  in  previous  shock 
tube  experiments   'L   ,  but  the  initial  conditions  to  determine  it, 
where  again  different,  in  this  case,  the  calibration  runs  were 
accomplished  in  the  range  from  1  to  10  atmospheres  of  pressure. 

As  a  preview  of  the  following  work,  Chapter  II  describes  the 
equipment,  instrumentation  and  procedures  used  in  the  course  of  this 
work.  Chapter  III  is  related  to  the  general  theory  behind  shock  tube 
operation  as  well  to  the  theory  related  to  the  experiment  accomplished. 
Chapter  IV  is  devoted  to  the  experiment  itself,  the  data  reduction 
and  presentation  of  results,  and  Chapter  V  gives  the  conclusions  of 
the  present  work. 

It  is  proper  to  hereby  acknowledge  the  kind  assistance  of  several 
people  who  contributed  to  the  completion  of  this  work.  The  author 
wishes  to  express  his  appreciation  to  his  advisor,  Dr.  Daniel  J.  Collins, 
for  his  guidance  and  assistance  during  the  entire  time  this  work  was 
carried  out;  to  Mr.  Dana  Maberry  and  Mr.  Norman  Leckenby  for  their 
cooperation  and  assistance  in  the  performance  of  the  experimental 
program;  to  Mr.  Robert  Besel  for  his  assistance  with  the  graphics  work. 

Last,  but  not  least,  the  author  is  particularly  indebted  to  his 

wife,  Maria  Angelica,  for  her  assistance  in  the  preparation  of  the 

manuscript  and  for  her  understanding  and  patience  during  the  entire 

course  of  studies. 
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CHAPTER  II 
DESCRIPTION  OF  EXPERIMENTAL  EQUIPMENT 
I.  SHOCK  TUBE 

General 

The  shock  tube  shown  in  Plate  I  was  constructed  at  the  Naval 
Postgraduate  School,  Monterey,  California.  The  equipment  has  been 
installed  in  one  of  the  Aeronautical  Laboratories  of  the  School. 

Shock  tube  configuration  was  chosen  with  two  definite  purposes 
in  mind,  namely  to  be  used  as  a  laboratory  equipment  for  the  courses 
in  the  Gas  Dynamics  option  in  the  Aeronautical  Engineering  curricula 
at  the  School  and  at  the  same  time  as  a  research  aid  in  the  field 
of  high  temperature  gas  dynamics. 

Essentially  a  shock  tube  consists  of  a  rigid  cylinder  (although 
different  cross-sectional  shapes  are  possible),  divided  into  two 
sections  by  a  thin  diaphragm  mounted  normal  to  the  axis  of  the  tube. 
The  principle  of  operation  calls  for  a  pressure  difference  to  exist 
across  the  diaphragm.  The  high-pressure  side  is  called  the  Driver 
section  and  a  low  pressure  region  is  called  the  Driven  section. 

Figure  I  shows  a  schematic  of  the  shock  tube  sections  concerned, 
The  material  used  is  a  stainless  steel  tube,  3  in.  in  internal 
diameter,  with  a  wall  thickness  of  0.5  inc.  and  a  highly  nolished 
internal  surface.  The  overall  length  of  the  tube  is  approximately 
25  ft.  constructed  in  four  sections  of  different  lengths.  This 
arrangement  facilitated  the  construction  of  the  tube  and  at  the  same 
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PLATE      I 

SHOCK     TUBE     AT     NAVAL      POSTGRADUATE      SCHOOL 

MONTEREY,    CALIFORNIA 
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time  provides  flexibility  in  the  case  that  different  driver  or  driven 
tube  lengths  are  required.  The  first  section  of  tubing,  with  a  length 
of  7  ft.  3  in.  is  the  driver  section.  The  other  three  pieces  of  the 
tube  make  up  the  driven  section.  Between  the  driver  section  and  the 
driven  section  is  inserted  a  transition  section  which  holds  the  valve 
body  and  diaphragm  section.  Connecting  the  three  tube  sections  that 
make  up  the  driven  end  are  two  6  in.  outside  diameter  instrumentation 
sections  of  4.5  in.  length.  These  have  two  instrumentation  ports 
at  each  side  of  the  tube  and  two  on  the  upper  part,  intended  to  house 
the  different  types  of  detectors  which  can  be  used,  such  as  pressure- 
transducer,  shock  detectors,  temperature-sensitive  resistance  (thin 
metal  film),  ionization  detectors,  light  screen  schlieren,  etc. 

The  end  of  the  shock  tube  has  another  instrumentation  station, 
with  four  ports  located  90  ° apart  of  each  other.  The  end  plate  of 
the  tube,  also  has  provision  for  the  insertion  of  instrumentation. 
All  instruments  and  plugs  used  in  the  instrumentation  sections  are 
flush  mounted  with  the  inside  wall. 

Different  sections  are  bolted  together  and  end  flanges  and  joints, 
are  sealed  by  means  of  Viton  "0"  rings.  The  ratio  between  the  length 
of  the  tube  and  the  inside  diameter  is  76  and  the  ratio  between  the 
length  of  the  driver  and  driven  section  is  0.368. 

As  indicated  before,  the  driver  and  driven  section  are  separated 
by  a  diaphragm.  In  this  shock  tube,  there  are  two  diaphragms  separated 
by  a  distance  of  2.5  in.  called  the  diaphragm  chamber,  which  will 
be  described  later,  permits  a  very  ingenious  but  simple  way  of  ob- 
taining a  rather  precisely  controlled  bursting  pressure.  The  shock 
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tube  assembly  is  supported  as  shown  in  Plate  I. 

Driver  End 

A  very  important  consideration  in  the  operation  of  the  shock 
tube  is  the  vacuum  and  gas  filling  system  required.  Figure  2  shows 
a  schematic  of  the  system  used  for  the  driver  end  section  of  this  tube. 
The  vacuum  pump  evacuates  the  high  pressure  region  of  any  remaining 
gas,  air  and  moisture  before  the  section  is  filled  with  the  gas  to  be 
used.  Gas  is  admitted  to  the  front  end  section  of  the  shock  tube  from 
a  high  pressure  bottle  of  gas  by  means  of  tubing  through  a  regulator 
valve  connected  directly  to  the  bottle  and  a  pressure  control  valve 
located  beside  the  tube.  By  means  of  a  "T"  connection  two  high  pressure 
bottles  of  either  the  same  or  different  gases,  can  be  used  in  the 
charging  of  the  shock  tube.  The  pressure  of  the  gas  in  the  tube  is 
obtained  from  a  pressure  gauge  located  in  the  upper  part  of  the  front 
end  section  of  the  tube.  To  provide  a  means  of  evacuating  the  gas, 
after  the  shock  tube  has  been  fired,  a  bleed  valve  is  located  in  the 
bottom  part  of  the  front  end  section.  Plate  2  shows  the  actual  set-up 
of  the  driver  end  section. 

Driven  End 

The  vacuum  system  of  low  pressure  region  or  driven  end  is  an  impor- 
tant part  of  shock  tube  operation.  Figure  3  shows  a  schematic  drawing 
of  the  different  parts  which  make  up  the  vacuum  system.  Its  details 
and  purposes  are  described  below.  Basically  the  main  concern  faced  in 
the  design  of  the  vacuum  system,  is  the  range  of  vacuum  required. 
With  this  criteria  in  mind,  the  selection  of  a  system  which  will  be 
adequate  and  expeditious,  is  of  prime  importance  for  optimum  utilization 
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of  the  shock  tube. 

Connected  directly  to  the  shock  tube  valve  body  is  the  main 
valve,  which  controls  the  operation  of  the  vacuum  system.  The  main 
valve  is  connected  to  a  General  Electric,  Model  5KC45PG4X,  vacuum 

pump  which  is  designed  to  initially  evacuate  the  driven  end  and  is 

_3 
capable  of  pulling  up  a  vacuum  of  approximately  10   mm.  of  mercury. 

Also  connected  to  the  main  valve  is  a  diffusion  pump,  which  is  canable 
of  pulling  up  vacuums  of  10"   torr.  The  diffusion  pump  used,  is 
manufactured  by  National  Research  Corporation  (NRC),  a  subsidiary  of 
Norton  Company  of  Newton,  Massachusetts.  The  system  also  consists  of 
a  Circular  Chevron  Cryo-Baffle  and  a  pneumatically  operated  High  Vacuum 
Slide  Valve,  both  manufactured  by  NRC.  The  equipment  is  shown  schemat- 
ically in  figure  3.  Two  thermocouple  vacuum  gauges  and  one  ionization 
gauge  tube,  are  used  to  measure  the  pressure.  The  thermocouples 
gauges  are  manufactured  by  NRC  and  the  type  used  is  0531.  The  vacuum 
ionization  control  gauge  instrument  is  also  manufactured  by  NRC  and 
its  model  is  NRC  831.  The  ionization  gauge  tube  is  a  RG  75-K  manu- 
factured by  Vacuum  Electronics  Corporation  (Veeco),  Long  Island,  New 
York.  The  vacuum  gauge  instrument,  vacuum  ionization  gauge  instrument 
and  the  control  switches  for  the  diffusion  pump,  vacuum  pump  and  pneu- 
matic actuator  for  the  high  vacuum  slide  valve  are  mounted  in  an  instru- 
ment panel  for  easier  control  and  operation  of  the  vacuum  system. 
Plate  3  shows  the  actual  set-up  of  the  driven  end  vacuum  system  and 
instrument  panel  respectively. 
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PLATE      3 
DRIVEN      END     VACUUM      SYSTEM 
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II.   INSTRUMENTATION 

General 

The  successful  use  of  the  shock  tube  depends  entirely  on  suit- 
able instrumentation  for  the  required  flow  measurements.  Although  the 
type  of  instruments  to  be  used  depend  largely  on  the  experiment  being 
run,  accuracy  needed,  and  money  available,  they  all  require  a  precise 
and  fast  response,  since  the  duration  of  the  shock  wave  in  the  tube 
is  of  the  order  of  milliseconds.  Consideration  is  given  here  to  the 
actual  instrumentation  used  in  the  shock  tube  at  the  Naval  Postgraduate 
School . 

Temperature  sensitive  resistance  detectors 

One  of  the  fundamental  measurements  in  shock  tube  tests  is  that 
of  the  shock  speed.  One  way  to  accomplish  this  measurement  is  by 
using  temperature  sensitive  resistance  detectors  or  most  commonly 
known  as  thin  film  gauges.  These  devices  are  of  local  manufacture, 
using  a  short  length  of  glass  or  pyrex  about  0.25  in.  in  diameter,  with 
one  flat  end  and  edges  carefully  polished.  A  film  of  Liquid  Bright 
Platinum  No.  05,  manufactured  by  Engelhard  Inc.,  East  Newark,  New 
Jersey  was  used  and  a  thin  strip  was  painted  over  the  flat  end  and 
extended  through  the  sides  of  the  tube;  the  film  should  be  around  1/32" 
in  width  and  its  length  and  thickness  such  that  the  resistance  is  about 
100  ohms;  after  that  is  accomplished  it  is  let  to  air  dry  for  a  minimum 
of  one  hour,  then  baked  in  an  oven  to  1240°F  for  20  minutes  and  then 
air  cooled;  the  process  is  completed  by  annealing  three  times  to  1100°F 
for  1  hour,  allowing  complete  air  cooling  every  time.  Leads  are 
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soldered  to  the  sides  of  the  gauge  by  using  Eccobond  Solder  V-91 ,  a 
product  of  Emerson  &  Cuming  Inc.,  Canton,  Mass.,  and  after  that,  the 
gauge  is  placed  on  a  phenolic  plug,  for  insulation  purposes,  to  hold 
it  in  place  in  the  instrumentation  port  plug  of  the  tube,  as  shown  in 
Figure  4.  Although  "0"  rings  are  provided  for  sealing,  in  order  to 
avoid  any  possibility  of  leak  through  the  thin  film  and  the  phenolic 
plug,  the  free  space  between  them  was  filled  up  with  Enoxy  made  with 
APCO  210  resin  and  APXO  180  hardneer,  products  of  Applied  Plastic  Co., 
Inc.  (APCO),  El  Segundo,  California. 

The  response  of  a  thin  film  gauge  operates,  dependent  on  a  change 
in  resistance  with  a  change  in  temperature  of  the  thin  metal  film. 
In  the  case  of  the  shock  tube,  this  rise  in  temperature  is  provided 
by  the  shock  wave.  With  a  constant  current  passed  through  the  metal 
film,  a  rapid  increase  in  voltage  is  obtained  as  the  shock  wave  passes 
over  the  gauge.  This  voltage  properly  amplified,  is  fed  into  an 
oscilloscope  and  recorded. 

Three  thin  film  gauges  were  used  for  the  measurement  of  the  shock 
wave  speed.  Each  one  of  them  was  mounted  in  one  of  the  instrumentation 
station  ports,  provided  for  that  purpose,  as  shown  in  Figure  1. 

Calibration  of  these  gauges  was  not  necessary  because  the  signal 
obtained  from  them,  was  used  to  trigger  the  oscilloscope  and  from 
their  stepwise  signal,  time  resolution  was  obtained  which  allowed 
determination  of  the  shock  speed. 

Heat  transfer  thin  film  gauges 

This  type  of  gauge  is  the  same  as  the  one  described  above.  A 
different  name  is  used  for  purpose  of  identification  since  its  intended 
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use  is  different.  The  principle  of  operation  is  the  same  as  the  one 
described  for  the  temperature-sensitive  resistance  detectors  and  the 
manufacture  process  was  done  in  the  same  manner.  This  gauge  is  located 
on  the  end  plate  of  the  Driven  end  of  the  shock  tube,  mounted  flush  with 
the  interior  surface  of  the  plate  and  held  in  place  by  means  of  a  phenolic 
plug  of  similar  construction  to  the  one  shown  in  Figure  4.  The  gauge 
measures  the  temperature  of  the  end  wall  before  and  after  the  arrival  of 
the  shock  wave.  The  signal  produced  by  the  gauge  is  fed  into  an  oscill- 
oscope and  is  recorded  photographically. 

It  was  necessary  to  calibrate  the  heat  transfer  gauge  in  order  to 
establish  the  temperature  coefficient  of  resistivity  (a).  The  cali- 
bration was  done  by  placing  the  heat  transfer  gauge  into  a  bath  of 
ice  water  and  recording  the  initial  gauge  temperature  and  resistance; 
the  water  was  then  slowly  heated  to  a  temperature  approximately  of 
90°  C,  and  readings  of  temperature  and  resistance  of  the  gauge  were 
recorded  during  the  process.  The  following  relationship1-  J,L  J  was 


used 


RF  -  R0  [1  +  a(Tp  -  T„)  ] 


where 

RF  =  gauge  resistance 

TF  =  gauge  temperature 

R   =  cold  gauge  resistance 

T   -  initial  temperature 


The  temperature  coefficient  of  resistivity  (a)  was  obtained  by 
plotting  the  change  in  resistance  versus  the  change  in  temperature. 
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The  resistance  of  the  thin  film  gauge  was  approximately  100  ohms, 
and  was  found  to  be  a  linear  function  of  its  temperature.  Several  runs 
were  made  for  the  calibration  and  one  example  is  given  in  Figure  5. 
The  plot  was  fitted  using  the  least  squares  method.  The  temperature 
coefficient  of  resistivity  was  found  to  be  0.00104  ohm/ohm  deg.  C. 

Electronics 

In  Figure  6  is  shown  a  schematic  of  the  electronic  set-up  used 
in  the  shock  tube.  Gauges  Nos.  1,  2,  and  3,  represent  the  thin  film 
gauges  used  for  the  measurement  of  temperature  difference  at  the  end 
wall.  The  gauges  are  connected  to  a  bridge  as  shown  in  Figure  7. 
Since  the  resistances  R-, ,  R«»  Rq  and  R*  are  much  larger  than  RF,  the 
current  I  through  the  gauge  is  maintained  steady,  thus  as  the  change 
in  gauge  resistance  is  attained  with  the  passing  of  the  shock  front, 
a  stepwise  signal  in  voltage  is  produced  and  this  signal  fed  into  the 
oscilloscopes.  In  the  arrangement  shown  in  Figure  6,  the  signal  from 
gauge  No.  1,  goes  to  an  Amplifier  where  it  is  amplified  about  two 
hundred  times  and  from  there,  is  used  to  trigger  both  oscilloscopes. 
The  output  from  gauges  No.  2  and  3  are  connected  to  a  dual  beam  oscill- 
oscope, where  the  distance  between  the  two  peaks  in  their  traces,  gives 
the  time  delay  needed  to  compute  the  shock  wave  speed.  The  output 
signal  from  the  H.T.  gauge  is  connected  to  the  second  single  beam 
oscilloscope.  The  electronic  equipment  used  is  a  Wide  band  Differ- 
ential D.C.  Amplifier  Model  885,  manufactured  by  Astrodata,  Inc., 
Anaheim,  California;  a  Type  551  Dual  Beam  Oscilloscope  and  a  Type  549 
Storage  Oscilloscope,  both  manufactured  by  Tektronic,  Inc.,  Portland, 
Oregon. 
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III.  OPERATING  PROCEDURES 

Filling  and  Vacuum  System 

The  first  step  in  operating  the  shock  tube,  is  to  obtain  a  vacuum 
in  the  Driven  end.  With  reference  to  Figure  3  and  starting  with  all 
valves  closed,  the  vacuum  pump  is  turned  on,  and  when  ready  to  pull 
down  vacuum,  control  valve  No.  2  and  the  Main  valve  are  open;  control 
valve  No.  1  connecting  the  vacuum  pumn  with  the  diffusion  pump,  is  a 
two-way  valve,  which  in  the  open  position,  closes  the  line  to  the 

diffusion  pump.  Working  in  that  condition  it  is  possible  to  obtain  a 

-3 
vacuum  of  approximately  10   mm.  of  mercury;  readings  are  obtained 

from  either  the  vacuum  gauge  and/or  thermocouple  gauges.  If  a  higher 

vacuum  is  needed,  the  diffusion  pump  is  used  and  one  opens  the  vacuum 

slide  valve  and  close  control  valves  No.  1  and  2.  In  this  condition 

the  system  is  operated  until  the  desired  vacuum  is  obtained.  Although 

it  was  not  tested  during  the  present  work,  with  the  diffusion  pump  it 

-9 
is  possible  to  obtain  vacuums  in  the  neighborhood  of  10  "  mm.  of  mercury, 

When  the  desired  vacuum  is  obtained,  the  main  valve  is  closed  and  then, 

by  means  of  the  quick  acting  filler  valve  the  working  gas  is  pumned 

into  the  Driven  end.  The  gas  is  usually  stored  in  high  pressure 

bottles. 

The  next  step  is  to  fill  the  Driver  end  with  the  gas  to  be  used, 

to  the  desired  pressure.  A  vacuum  pump  could  be  used  in  the  driven 

end,  in  order  to  assure  a  clean  and  moisture  free  atmosphere  before  the 

working  gas  is  admitted.  With  reference  to  Figure  2,  the  gas,  usually 

stored  in  high  pressure  bottles,  is  fed  into  the  driver  section,  by 
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means  of  the  pressure  control  valve  and  at  the  same  time  to  the  dia- 
phragm chamber.  Assuming  a  pressure  P*  is  desired  in  the  driver  end, 
the  diaphragm  chamber  is  filled  up  to  a  pressure  P\   less  than  P,. 
The  chamber,  actually  a  very   small  compartment  of  the  tube  in  length, 
is  sealed  off  by  two  diaphragms,  the  first  being  thick  enough  to 
withstand  the  pressure  difference  P-  -  Pi  and  the  second  one  to  with- 
stand the  pressure  Pi  but  not  P*.  When  the  required  pressure  has 
been  reached  in  the  driver  section,  the  diaDhragm  chamber  is  either 
rapidly  evacuated  or  expanded  and  the  two  diaphragms  are  then  burst 
almost  at  the  same  time.  This  rather  simple  method  of  runturing  the 
diaphragm  has  the  main  advantage  of  giving  a  precise  control  pressure 
at  which  the  diaphragms  are  ruptured  and  permits  firing  of  the  tube 
at  the  desired  moment.  This  ability  to  control  the  firing  of  the 
tube  is  highly  desired  if  for  any  reason  something  went  wrong  with 
the  electronics  and  data  recording  equipment  at  the  last  moment. 

To  evacuate  the  shock  tube  after  it  has  been  fired,  bleed  valves 
are  provided  along  the  sections  of  the  tube  as  indicated  in  Figures 
2  and  3. 

Data  recording 

The  successful  run  of  the  experiments  depends  entirely  in  the 
ability  of  recording  the  required  data.  Due  to  the  short  time  every 
shot  lasts,  the  low  output  signals  of  the  thin  film  gauges  and  the 
high  sensitivity  of  the  equipment  to  be  used,  this  part  of  the  experi- 
ment required  some  coordination  of  the  many  factors  involved. 

The  data  recorded  was  the  following:  driver  end  pressure, 
driven  end  pressure,  initial  temperature  in  the  shock  tube,  time 
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required  for  shock  waves  to  travel  a  fixed  distance  and  voltage  rise 
in  the  end  wall  at  shock  arrival.  The  first  three  measurements  were 
taken  by  direct  reading  of  the  instruments  provided  for  that  purnose. 
The  time  required  for  the  shock  wave  to  travel  a  fixed  distance  in 
the  shock  tube,  namely  the  distance  between  thin  film  gauges  No.  2 
and  3*-  1^ur   J  was  obtained  by  means  of  a  dual  beam  oscilloscope  and 
recorded  using  a  Polaroid  film.  A  typical  trace  of  the  signal  obtained 

is  shown  in  Plate  4.  With  above  data,  the  speed  of  the  shock  wave  was 

T2l 
determined  and  knowing  the  speed  of  sound  in  the  gas  usedL  J  the  shock 

wave  Mach  number  was  computed. 

The  temperature  rise  in  the  end  plate  is  given  by  the  following 

relationship: 


AT  =  T  -  T  =       aE 

W     0 


a  I  I  Z    P\ 
x   0  x   0 


The  temperature  coefficient  of  resistivity  (a)  was  previously 
determined  for  the  gauge  as  described  before;  the  current  (I  )  and 
resistance  (R  )  were  measured  before  the  runs.  The  change  in  voltage 
was  recorded  from  the  oscilloscope  hook  to  the  heat  transfer  thin 
film  gauge  shown  in  Figure  6.  A  typical  trace  of  the  signal  obtained 
is  given  in  Plate  5. 
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PLATE      4 

TRACE      OF     TEMPERATURE      SENSITIVE 
RESISTANCE      DETECTORS 


34 


PLATE       5 
TRACE     OF     HEAT      TRANSFER      THIN     FILM      GAUGES 
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CHAPTER  III 
THEORY  BEHIND  SHOCK  TUBE  OPERATION 

I.  GENERAL  THEORY 

The  theory  of  the  shock  tube  is  well  known  and  clearly  described 
in  many  textbooks.  The  same  comment  can  be  made  in  regard  to  the  theory 
of  shock  waves.  It  is  not  the  intended  purpose  of  this  thesis  to  go 
deeply  into  the  deviation  of  the  different  relations  governing  the 
actual  phenomenon  occuring  in  the  tube,  but  to  give  a" simple  presentation 
on  the  production  of  a  shock  wave  and  then  look  in  more  detail  at  the 
theory  governing  the  experiment  which  was  performed. 

In  Figure  8  is  shown  the  thermodynamic  history  of  the  shock  tube. 
Figure  8a  and  8b  corresponds  to  time  t  =  0.  The  Driver  end  corresponds 
to  region  (4)  and  the  Driven  end  to  region  (1).  At  t  =0  the  temperature 
of  both  regions  is  the  same. 

At  the  appropriate  time  the  diaphragm  is  broken.  A  shock  wave  (S) 
is  propagated  into  the  driven  end  and  a  rarefraction  (or  expansion)  wave 
(R)  is  propagated  back  into  the  driver  end.  This  is  shown  in  Figure  8c 
at  time  t-.  >  0.  As  the  shock  front  travels  into  region  (1),  the  fluid 
is  compressed  to  a  higher  pressure  P2  (Figure  8d)  and  higher  temperature 
T?  (Figure  8e)  in  a  nonisentropic  process.  On  the  other  hand,  as  the 
rarefraction  wave  travels  into  region  (4),  the  fluid  is  expanded  to  a 
lower  pressure  P^  and  a  lower  temperature  T~  by  a  process  which  is 
essentially  isentropic.  Thus,  although  the  pressures  Pp  and  P-  are  equal, 
there  are  two  regions  between  (R)  and  (S).  These  two  regions  correspond 
to  regions  (2)  and  (3)  in  the  graphs  of  Figure  8.  The  gas  in  region  (2) 
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has  been  compressed  and  heated  by  the  shock  wave,  however,  the  gas  in 
region  (3)  has  been  expanded  and  cooled.  The  narrow  zone  called  the 
contact  surface  or  contact  discontinuity  indicated  by  (C)  in  Figure  8c 
represents  the  line  where  a  temperature  and  density  discontinuity  exists, 
although  pressures  are  the  same  on  each  side  of  (C).  Physically  it 
denotes  the  position  occupied  by  that  gas  which  was  originally  at  the 
diaphragm  and  is  the  point  where  the  experimental  gas  and  the  driver  gas 
make  contact.  This  "contact  surface"  travels  along  the  tube  behind  the 
shock  front.  The  above  explanation  is  the  ideal  one.  In  the  real  case 
some  mixing  of  gases  is  usually  present  at  the  contact  surface  which 
causes  the  temperature  drop  to  occur  in  a  finite  time  period. 

The  most  common  notation  for  the  different  regions  is  the  one  used 
in  the  above  explanation.  Another  region,  called  region  (5),  is  present 
in  the  shock  tube  operation  where  the  shock  wave  is  reflected  at  the  end 
wall  of  the  tube  with  a  further  rise  in  temperature  and  pressure  exists. 
This  region  is  shown  in  Figure  9,  which  actually  represents  a  physical 
plot  (distance  vs  time)  of  the  phenomenon  initially  produced  in  the  shock 
tube  after  the  diaphragm  has  been  burst.  It  can  be  seen  in  this  plot, 
that  the  rarefraction  (or  expansion  wave)  covers  a  zone  characterized  by 
falling  temperature,  pressure  and  density  below  the  initial  conditions. 
These  waves  spread  with  time  and  the  fluid  properties  changes  occur  in  a 
more  gradual  process  as  compared  with  the  shock  front  and  contact 
surface.  Due  to  their  spreading  action,  these  waves  are  called  the 
expansion  fan,  with  the  initial  and  terminal  wave  called  the  head  and 
tail  waves  respectively. 
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The  theoretical  analysis  for  shock  tube  phenomenon  is  greatly 
simplified  by  the  assumption  of  one-dimensional  flow,  where  the  flow 
properties  are  considered  to  be  uniform  over  any  cross-section  and  to 
vary  only  with  distance  along  the  tube  and  in  general  with  time. 

The  fundamental  equations  governing  a  non-stationary  one-dimensional 
flow  are  the  continuity,  momentum,  energy  and  equation  of  state  and  they 
are  given  by  the  following  relationships  respectively: 


(1) 
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The  most  common  procedure  to  solve  for  the  relationships  aoplied  to 

shock  tube  operation  is  using  the  above  equations  together  with  the 

methods  of  characteristics.  A  detailed  coverage  of  these  derivations  is 

given  in  Reference  [1]  and  [2],  and  for  the  method  of  characteristics 
itself  in  Reference  [9]. 
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Equation  (5)  relates  the  initial  applied  pressure  ratio  across  the 
diaphragm  to  the  shock  wave  Mach  number.  It  gives  a  measureing  of  the 
strength  of  the  shock  in  a  given  test  gas  and  shows  how  for  a  given 
pressure  ratio,  the  shock  strength  depends  on  the  specific-heat  ratio 
and  sound  speed  in  the  high  pressure  gas. 

The  strongest  shocks  can  be  obtained  by  using  a  driver  gas  having 
a  high  speed  of  sound  and  low  specific  heat  ratio. 

Equations  (6)  and  (7)  are  useful  to  determine  the  state  of  the  gas 
behind  the  reflected  shock.  In  the  form  presented  they  are  specially 
suited  for  the  experiment,  leading  to  a  full  solution  of  the  pressure 
and  temperature  behind  the  reflected  shock  in  terms  of  the  incident  Mach 
number  and  initial  pressure  and  temperature,  parameters  measured  directly 
from  the  experiment. 

II.  THEORY  FOR  EXPERIMENTS 
The  region  behind  the  reflected  shock  wave  (region  5)  is  idealized 
to  consist  of  a  hot  semi-infinite  gas  adjacent  to  a  semi-infinite  solid. 
Viscouse  dissipation  may  be  neglected  and  the  pressure  assumed  constant 
in  the  gas  boundary  layer  near  the  end  wall  of  the  tube1-  .  Under  these 
conditions,  we  may  write  the  equations  of  continuity  and  energy  as 


(«)      -it    +  -53T  <*>">  *  ° 


and 


(9)  p<=p<-ir  +  u^>  -  ~i^^r) 
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The  continuity  equation  is  satisfied  by  introducing  a  stream 
function  \p,   such  that: 


(10) 


j) _   3^  pu   _     d\j) 

p       dx  p         8t 

Mw  Hw 


Changing  coordinates  from  x  to  tp,  the  energy  equation  becomes 


(id        -^-(-H-)    ■  -sS- ( 481  x-IS- ) 


a,  3t     '  8^     v     9  8ip 

w  ip 


where 


(12)  %    -•  J&*       ,  e  -    -f-     ,  K(e)    =   -f    ,  -e_   =    -ft 

w  pwcp  w  kw  pw  T 


With  the  above  transformation,  no  characteristic  time  appears  in 

2 

the  equation,  leaving  6  as  a  function  of  the  expression  — *j —  ,  implying 

V 

that  the  wall  surface  temperature,  T  ,  is  a  constant. 

w 

Introducing  a  similarity  parameter 


(13)      n  = 


<  V> 1/2 


The  one  dimensional  time  dependent  energy  and  continuity  equations  near 
the  end  wall  of  the  tube  reduce  to 


»«      Tf<^*-t->  +  ^  =  0 


subject  to  the  following  boundary  conditions: 

T 

(15)     e  (0)  =1.0   ,  e  («)  =  — 

T 
w 
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The  subscripts  w  and  °°  refer  to  numerical  values  evaluated  at  the 
end  wall  condition  just  after  the  shock  wave  arrival  and  to  any  quantity 
measured  behind  the  reflected  shock  wave,  respectively. 

Equation  (14)  can  be  integrated  from  n  =  0  to  n  =  °°  ,  after  the 
temperature  dependence  of  the  thermal  conductivity  is  established. 
Reference  [5]  and  [6]  assume  the  thermal  conductivity  to  vary  as  a  power 
law  in  temperature 

With  that  assumption,  equations  (14)  and  (15)  can  be  used  to 
obtain  the  following  relations: 

TO     e-  "  ^   -    «*%)      •     %  -   (-f-)  n=0 

By  equating  the  heat  flux  at  x  =  0  to  the  heat  flux  at  the  surface 

li&l 
of  the  solid  end  wall  it  can  be  shownL  J  that 


«     %-W^ 


where 


(19)       A  = 


1/2 


ks  ps  cs 
ew    I  kw  pw  cp 


In  the  above  equations,  the  subscripts  s  and  w  refer  to  properties 
of  the  gauge  and  properties  of  the  gas  at  the  wall  temperature  respectively, 
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Since  the  quantities  0^  and  q  can  be  experimentally  evaluated 
from  each  run,  a  series  of  tests  yields  an  experimental  relation  between 
them  for  the  gas  tested.  Equation  (16)  established  that  the  relation  of 
these  quantities  is  dependent  on  the  exponent  a.  Then,  by  comparing  the 
theoretical  and  experimental  data  it  is  possible  to  choose  a  value  of 
"a"  which  best  fits  the  data.  The  numerical  value  of  "a"  was  computed 
by  minimizing  the  sum  of  the  squares  of  the  deviation  between  the 
experimental  and  theoretical  curves.  Once  it  was  established,  the 
temperature  dependence  of  the  thermal  conductivity  is  given  by  equation 
(16). 
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IV_   EXPERIMENTAL  RESULTS  AND  DISCUSSION 
I.   DESCRIPTION  OF  EXPERIMENT 

Purpose 

The  purpose  of  the  experiment  described  herein  was  to  look  at 
thethermal  conductivity  of  argon  at  relatively  higher  pressure,  i.e. 
10  to  50  atmospheres  and  to  compare  the  results  obtained  with  those 
similar  experiments,  done  at/or  around  atmospheric  pressure.  Previous 
data  for  experiments  having  the  same  conditions  of  the  one  done  for 
this  thesis  was  not  available,  so  a  comparison  in  that  sense  was  not 
possible. 

Calibration  Runs 

With  reference  to  the  theory  and  equations  presented  in  the  pre- 
vious chapter,  it  was  necessary  to  calibrate  the  thin  film  heat  transfer 
gauge  to  obtain  a  value  of  B_>  in  order  to  compute  later  the  experi- 
mental values  of  q  .  This  calibration  was  done  in  the  shock  tube, 
w 

by  making  a  series  of  runs  (30  in  total),  where  argon  was  heated  to 
temperatures  ranging  from  1000  to  1500  degrees  Kelvin.  In  this  range 
of  temperature,  the  thermal  conductivity  of  argon  has  been  experi- 
mentally determined.  To  determine  the  heat  transfer  rate  q  ,  a  value 

J  w 

of  the  parameter  "a"  is  required.  The  value  a  =  0.697  taken  from 

Collins  and  Menard*-  *   was  used.  Equation  18  was  then  used  to  obtain  a 

?  1 1? 

value  for  6  =  0.062060  cal/cm  deg.  K  (sec)  '  with  a  standard 

deviation  of  +  0.010.  This  value  was  found  to  be  around  15%  higher 

than  others  previously  determined.  Due  to  this  reason,  and  that  the 
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calibration  runs  were  accomplished  also  at  relatively  high  pressures, 
between  1  and  10  atmospheres,  the  validity  of  the  use  of  the  parameter 
a  =  0.697  is  questionable.  Unfortunately,  due  to  lack  of  time  it 
was  not  possible  to  go  into  the  determination  of  the  optimum  value  of 
"a"  to  be  used  in  this  case. 

Experimental  Runs 

The  experimental  runs  were  made  in  the  shock  tube  previously 
described  using  the  data  recording  equipment  and  instruments  mentioned 
before.  The  gas  tested  was  Argon  using  Helium  as  a  driver.  The  runs 
were  made,  as  to  obtained  temperatures  between  1500° to  5000° K  and 
pressures  of  10  to  50  atmospheres  in  the  region  behind  the  reflected 

shock  wave  (region  5).  In  the  determination  of  thermal  conductivity, 

r  on 

the  value  of  k  was  obtained  from  the  values  given  in  CookL  .  Thus 
w  3 

the  problem  was  reduced  to  find  a  value  of  "a"  for  equation  (20)  to 
fit  the  data. 

(20)    K  (e)  =   k    =  oa 

w 

A  list  of  data  recorded,  as  well  as  the  values  computed  is  given  in 
Plate  6. 

II.  DATA  REDUCTION 
For  every  run  made,  the  data  obtained  was  the  driven  end  pressure 
(p-,),  the  shock  tube  temperature  (T  ),  the  time  resolution  for  the 
shock  wave  mach  number  computation  and  the  voltage  rise  at  the  end 
wall  thin  film  heat  transfer  gauge,  used  for  computations  of  temperature 
rise  in  the  end  wall.  With  that  data,  it  was  then  possible  to  obtain  a 
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dimensionless  plot  of  0  =  T  /  T  vs  (de/dn),,  =  q,  for  the  experiment 

r         00     00     yi  W     W 

accomplished.  After  that  was  done  and  by  using  equation  (14)  chapter 
III  and  by  comparing  the  theoretical  with  the  experimental  values 
obtained  and  minimizing  the  sum  of  their  deviations  a  plot  of  the 
residuals  versus  "a"  was  obtained  and  from  it,  the  best  value  of  "a" 
to  fit  the  data  was  selected. 

In  appendix  I  are  given  sample  computations  indicating  the 
breakdown  for  equation  (14)  obtained  in  order  to  solve  it  using  the 
IBM  360  computer  of  the  school  and  for  the  determination  of  q 
for  the  experiments  run. 

Although  several  partial  computer  programs  were  written, 
Appendix  II  shows  the  final  computer  program  written  that  solves 
for  the  best  value  of  a,  having  as  inputs  the  initial  data  obtained. 
The  output  of  those  programs  was  divided  in  three  steps;  the  first,  a 
tabular  array  given  for  each  run  the  initial  data  and  the  computed 
values  obtained  for  the  mach  number,  temperature  raise  in  the  end  wall, 
pressure  and  temperature  in  the  region  behind  the  reflected  shock, 
final  temperature  in  the  end  wall,  a  value  of  0  =  T/T  and  a  value  for 
q  .  With  these  last  two  quatnities,  a  subroutine  called  W4PL0T  was 
used  to  obtain  a  graphical  representation  of  the  experimental  results. 
Then,  equation  (14)  was  solved,  using  a  subroutine  called  SYSZDE  and 

from  it,  a  theoretical  value  of  0  =  T  /T  was  obtained,  for  different 

00    00'  yi  » 

values  of  "a".  The  residuals  between  the  theory  and  experiment  were 
then  obtained  and  these  values,  plotted  versus  a,  yielded  the  selection 
of  the  optimum  value  of  "a".  This  was  represented  in  a  granh,  using 
again  the  subroutine  W4PL0T.   The  final  plot  obtained  is  shown  in 
Figure  10. 
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PLOT    OF    EXPONENT   a  VS.    RESIDUAL 
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With  the  computer  program  shown  in  appendix  III,  a  graph  of  the 
experimental  data,  together  with  the  curve  for  best  fit  was  then  ob- 
tained and  at  the  same  time,  theoretical  curves  for  other  values  of 
"a"  chosen  arbitrarily  were  plotted  for  comparison,  as  shown  in  Figure 
11. 

Due  to  the  relatively  low  time  resolution  used  in  the  oscillo- 
scope connected  to  the  heat  transfer  thin  film  gauge,  1  millisecond 
vs  10  microseconds  of  reference  6,  the  flat  plateau  at  the  beginning 
of  the  trace,  as  the  one  mentioned  in  that  reference,  was  not  visible 
in  all  the  pictures. 

In  approximately  one  third  of  the  data  a  flat  plateau  was  iden- 
tifiable. Another  reference  point  visible  on  all  the  traces  was  used. 
This  plateau  point  was  some  25%  higher  than  the  reference  point.  A 
correction  factor  of  1.25  was  then  used  for  all  the  data. 
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CHAPTER  IV 

CONCLUSIONS 
The  experimental  data  is  preliminary  in  nature,  nevertheless 
the  thermal  conductivity  of  Argon  at  higher  pressures  10  to  50  atmospheres 
has  been  determined.  In  the  relationship  given  by  equation  16,  the 

experimental  value  of  a  =  0.4  was  obtained  with  k  computed  as  shown 

r  w 

in  appendix  II,  step  AKW.  No  previous  data  in  this  range  of  pressures 
has  been  obtained  so  that  it  is  not  possible  to  make  comparisons. 

It  is  necessary  perhaps,  in  some  future  study  of  the  thermal 
conductivity  of  Argon,  to  explain  from  a  theoretical  view  point  the 
decrease  in  the  value  of  "a"  with  pressure,  from  that  obtained  at 
lower  pressures. 

It  would  appear  that  it  is  also  necessary  to  repeat  some  aspects 
of  the  experiment  in  order  to  verify  some  of  the  results  obtained. 
This  applies  particularly  to  the  value  of  6  obtained. 

The  purpose  of  this  thesis  was  to  develop  a  fully  functioning 
shock  tube  and  to  demonstrate  the  functioning  of  the  shock  tube 
by  means  of  an  experimental  problem. 

The  shock  tube  is  now  a  working  tool  in  the  Aeronautical  Engin- 
eering Department  and  some  preliminary  results  have  been  obtained 
in  an  experiment  to  determine  the  thermal  conductivity  of  Argon  in 
the  pressure  range  from  10  to  50  atmospheres. 
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APPENDIX  I 
SAMPLE  CALCULATIONS 
a)  Solution  of  equation  (14)  for  data  reduction 


04)      tH^-SH  +  ^=  0 


from  equation  (12)  Chapter  III  is  found: 


K(9)  =  ea 


Replacing  in  equation  (14): 


2 
,a-l  d  6       d  /Qa-K  de       de      n 

.2      dn        dri       dn 


reducing  above  expression: 


/9,x       d20    a-1   de  2     Qa-1  de 
(21)       ^2  +  —   6T   +n  6   ^ 


Equation  (21)  represents  a  second  order  nonlinear  differential  equation 
In  order  to  solve  it,  it  was  reduced  to  a  system  of  two  first  order 
differential  equations  by  the  following  transformation: 

Let         X]  =  6     ,  x2  =  -[ji- 

Then,  the  final   equations  are: 
(22)  X.,'  =     x2 


(23)  x2'   =     1^    x22  -  nx/-9  x2 
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b)  Determination  of  q  for  experiment  runs 


The  value  of  3  was  known  to  be 


cal . 


3   =  0.06026 


cm2  K°(sec)1/2  J 


From  equations  (17),  (18)  and  (19)  Chapter  II,  is  obtained 


_de_ 

dn 


h      Tw"Ts 


IT    3, 


w 


w 


let 


AT  =  T  -  T 
w   s 


Using  equation  of  state  for  a  perfect  gas 


^./orr    tj ^^s. .  h  v  i  a 


w  w 


WW  w  ^w  p 


R  T        \  I   W   W   W   n 
W 


For  above  expression,  the  following  holds 


P   =  P    in  units  of  [cal /cm  ] 
w     °° 

cal 
k   is  given  in  units  of  [ 

w  cm°K  sec 

T   is  given  in  [°K] 


] 


P   =  (P  mm  Hg)  (3.184  10"15  ^fjg  )  is  a  conversion  factor 

c 

■^-    =     2.5  for  Argon 
R 


Then 


V 


de 
dn 


-=-  x  AT  x  3 

IT  S 


2.5  x  k  x  T  x  P  ' 

w    w    °° 
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